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Abstract. We present a review of the pseudo-SU(3) shell model and its application to heavy 
deformed nuclei. The model have been applied to describe the low energy spectra, B(E2) and 
B(M1) values. A systematic study of each part of the interaction within the Hamiltonian was 
carried out. The study leads us to a consistent method of choosing the parameters in the model. 
A systematic application of the model for a sequence of rare earth nuclei demonstrates that an 
overarching symmetry can be used to predict the onset of deformation as manifested through 
low-lying collective bands. The scheme utilizes an overarching sp(4,R) algebraic framework. 



1. Introduction 

Over the past years many calculations were carried out by using the SU(3) shell model and the 
pseudo-SU(3) symmetry. Ongoing improvements and innovations in experimental techniques (4ir 
detectors, radioactive beams, etc.) anticipate the identification of additional new phenomena 
in the near future, and more information about the existing ones. However, it is commonly 
accepted that the nuclear shell model should be able to address most of these issues and provide 
answers. The problem is that most shell-model theories are limited by the large dimensions 
of the required model spaces, which increases unmanageably with the number of nucleons A. 
Furthermore, more microscopic calculations are needed exactly in the region of heavy nuclei, in 
the mass range A > 100. 

Having this problem, it is essential to take full advantage of symmetries, those that are exact, 
as well as those which are only partially fulfilled. The selection rules that are associated with 
such symmetries generate, respectively, weakly coupled and disconnected subspaces of the full 
space and this allows for a significant reduction in the dimensionality of the model space. The 
symmetry used here is the pseudo-spin symmetry. 

A shell model theory for heavy nuclei requires a severe truncation of the model space. To 
reproduce the essential physics found in the low-energy states of a large space in a smaller one, 
we have to select the basis states relative to those parts of the interaction that dominate the 
low-energy structure. Nuclear physics supports the view that the nuclear effective interaction 
appropriate to low-energy excitations must have a strong correlation with the pairing and 
quadrupole-quadrupole interactions. The Q ■ Q interaction, which dominates for near mid- 
shell nuclei, is known to introduce deformation and led to the introduction of the SU(3) shell 
model. This led also to a very natural way to truncate large model spaces, namely to consider 
few basis states in correspondence with the eigenvalue of Q ■ Q operator. 



Based on the early contribution of Bohr and Mottelson[lj, Elliott[2] developed an SU(3) 
symmetry truncation for light ds-shell nuclei. Raju, Draayer, and Hecht[3j developed a pseudo- 
SU(3) scheme for heavy nuclei. Ground and gamma band rotational structure of deformed nuclei 
were described in 1984 by Draayer and Weeks [1], in the framework of the pseudo-SU(3) model 
by using only one SU(3) irreducible representation (irrep), the one corresponding to the highest 
deformation, and considering the abnormal parity spaces to seniority zero configurations. The 
fits to the available data were very accurate (\E exp — Ef^y < 25keV). In this small space the 
ground and the gamma (K=2) bands were very well described, as well as the inter-band and 
intra-band B(E2) strengths. The model used the mathematical and computational advances 
made by Draayer and Williams[5j, Akiyama and Draayer[6], Draayer, Weeks, and Hecht[7j, and 
Draayer and RosensteelJB]. 

Using mathematical and computational developments to the pseudo-SU(3) model, Castanos, 
Draayer, and Leschber improved it further to accommodate an increased model space, taking 
into consideration several protons and neutrons SU(3) irreps, as well as the total coupled 
representation in 1987[9l [TO] . In this version, many bands can be calculated and as well the 
B(E2) and B(M1) transitions values within and between bands. Strong B(M1) transitions were 
predicted first in various Gadolinium, Erbium, Uranium, and Plutonium isotopes using the 
pseudo-SU(3) model in 1987[9j. Naqvi and Draayer introduced a shell- model operator for K- 
band splitting in [11]. The K 2 operator is used to resolve multiple occurrences of L- values 
in SU(3) representations. In 1994, Bahri and Draayer completed an improved code for faster 
calculations [12]. Later the effects of paring interaction on the K-band mixing were studied[l~3Tll4j 
by Trotelnier, Bahri, Escher and Draayer within the framework of the pseudo-SU(3) model. 

Considering such larger spaces with an SU(3) Hamiltonian plus proton and neutron single 
particle energies, and pairing interactions, an extended version of the model was applied [T5l [To] 
to describe the low-lying energy structure and electromagnetic transitions probabilities of even- 
even nuclei in the rare-earth region . 

The purpose of this work is to show that by using a realistic Hamiltonian based on proton and 
neutron single-particle energies, Q ■ Q, and pairing interactions we can find systematically a set 
of parameters [16J, in the chosen model space, to obtain a good agreement with the experiment 
in the even-even heavy nuclei. We are also studying the effect of each term in the Hamiltonian 
plus the rotor terms, namely J 2 , K 2 , and Q3, on the energy spectra. The parameters of the 
rotor terms are fixed as the result of this study. 

The role of these parameters is further revealed in a systematic study of the behavior of first 
excited K n = 0;j~, and K 71 = 2 + bands in deformed even-even nuclei of the rare-earth region. 
The nuclei are ordered in F-spin multiplets of a Sp(4,R) classification scheme [17j. The energy 
levels of the ground band (g.b.) J = 2f , first excited J n = OF , , J n = 21, + states of 

nuclei that belong to the Fq = = 1/2(N 7T — N v ) multiplet (where N n and N v are the number 
of valence proton and neutron pairs) are plotted in Fig. [3j The energies of the same set of 
levels for nuclei belonging to the Fq = 1 multiplet are plotted in Fig. [4] Nuclei which belong to 
the Fq = multiplet have equal numbers of valence proton and neutron pairs, and those which 
belong to the Fq = 1 multiplet vary by two pairs of protons. With such an investigation [18], 
we provide a microscopic interpretation of the complex and varying behavior of these collective 
excitations, by applying the algebraic shell model with pseudo SU(3) symmetry. 

In the next section, the structure of the basis states is defined. An example is given for 156 Gd 
nucleus. The Hamiltonian is described in detail in section 3. The results of the analysis are 
explained in section 4. Also further applications of the model are presented. 



2. Basis states 

Building the many-body basis of the pseudo-SU(3) model relies on the proton-neutron degrees 
of freedom, pseudo-spin symmetry and group theoretical methods. The construction of the 
basis states starts with the deformed Nilsson single-particle levels for protons and neutrons, 
which are filled from the lowest energy upward, for a specific deformation parameter f3 for the 
corresponding nucleus. The completely filled shells are frozen and the remaining levels are taken 
into consideration. For heavy nuclei, the valence protons and neutrons fill different major shells. 
For a given nucleus there are two open shells (one for protons, one for neutron), each of them 
comprised of a set of normal and unique parity levels. The normal parity space is partitioned 
into pseudo SU(3) irreps, and the unique parity space is spanned by configurations of identical 
particles in a single j shell. It was shown that the seniority coupled scheme is appropriate 
for a description of the unique parity configuration [7]. Thus the low-energy structure of the 
normal parity part of the space will be dominated by a few irreps of pseudo SU(3). In the 
abnormal(unique) parity parts of the proton and neutron shells was shown that the configuration 
with low seniority are the most important ones, since those with high seniority are not favored 
(the pairing gap is large compared to the spacing of low- lying rotational bands.) 

In this way, the occupation numbers of the normal (N) and unique (C/)(also called intruder) 
parity states are determined for protons (it) and neutrons (u), , n%, , . For example, 
for 156 Gd, = 8 in % = 3, and n„ = 6 in % = 4 with degeneracies O,^ = 20 and U v = 30. 
Here, 7?7r[^], represent the pseudo associated numbers. For 8 protons in fj^ = 3 and 6 neutrons 
in fjv = 4 there are 32 possible SU(3) irreps in the proton space and 48 possible irreps in the 
neutron space, for proton and neutron pseudo-spin equal to zero. The corresponding irreps are 
ordered in descending C2 value and the first 6 proton and 5 neutron irreps are chosen, that 
are given in Table [T] The direct products of these proton and neutron irreps give a large set of 



Table 1. The SU(3) irreps for protons and neutrons corresponding to the largest C2 in Gd. 



(Att, ^w) 


c 2 


(K, fj>u) 


c 2 


(10,4) 


198 


(18,0) 


378 


(12,0) 


180 


(15,3) 


333 


(0,12) 


180 


(12,6) 


306 


(8,5) 


168 


(13,4) 


288 


(5,8) 


168 


(14,2) 


276 


(9,3) 


153 







possible strong coupled irreps from which we have selected the first 30 with the largest C% values. 
These irreps are given in Table [2] The many-particle states of n a nucleons in a shell (a = tt or u) 
belong to the totally antisymmetric irreps {1™"} of the unitary group of O,^ = (rj a + l)(r] a + 2) 
dimension. A complete classification of the basis states uses the quantum numbers following 
from subgroups of the unitary group. 

{fa} {fa} la (A«, Ma) S a K a L a J% 

U(n%) D U(tt% j2) x U(2) d SU(3) x SU(2) D 50(3) x SU(2) D SUj(2), (1) 

Above each group are given the quantum numbers that characterize its irreps with K a and 7 Q 
the multiplicity label of the indicated reduction. The eigenstates > of the Hamiltonian ([2]) 
are written as linear combination of the strong-coupled basis: 



(2) 



\<f>i >= |{ri 7r [/ 7r ](A 7r , /%), S^;n v [f v ]{X v , fj> v ), S u } p(X, p)KL S;JM >, 



(3) 



with two multiplicity labels: p that counts the number of times the irrep (A, p) occurs in the 
direct product (X^, p n )(X u , p u ) and K classifies the different occurrences of the orbital angular 
momentum L in (A, pi). 

Table 2. The first 30 pseudo SU(3) irreps used in the description of 15G Gd bands. 



(Att, Pit) 








total (A, p) 
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(25,7) 
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(27,3) 


(28,1) (23,8) (24,6) 


(10,4) 


(12,6) 


(22,10) 


(23,8) 


(24,6) 
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(10,4) 


(13,4) 


(23,8) 


(24,6) 
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(14,2) 


(24,6) 








(12,0) 


(18,0) 


(30,0) 


(28,1) 






(12,0) 


(15,3) 


(27,3) 








(12,0) 


(12,6) 


(24,6) 








(8,5) 


(18,0) 


(26,5) 


(24,6) 






(8,5) 


(15,3) 


(23,8) 


(24,6) 






(5,8) 


(18,0) 


(23,8) 









3. Model Interaction 

We studied in detail the effect of each term of the Hamiltonian on the energy spectra. The 
relative strength of each term in the interaction affects not only the energy spectra, but also the 
strength of the B(E2) and B(M1) transitions among those states. The Hamiltonian consists of 
quadrupole-quadrupole (Q -Q), proton (it) and neutron ( i/) single particle energies [H^ -\- Hgp) 7 
proton and neutron pairing ( Hp and Hp), as well as the rotor part (Hrot)- 

H = xQ ■ Q + + Hg p — G n Hp — G U H P + Hrot (4) 

where 

Hrot = aJ 2 + bK] + c 3 C 3 + a sym C 2 (5) 

C2 and C3 terms are the second and third order Casimir invariants of SU(3), which are related 
to the deformation of the nucleus. J 2 is the total angular momentum operator, and K 2 is the 
K-band splitting operator (resolving the multiple occurrences of the same angular momentum 
within an SU(3) irrep). The single particle energies are calculated in the standard way: 

H a sv = H + ds + dP = Hose - hO n(2Ls + pP) (6) 

where a = tt or v, Hq is the spherical oscillator for which SU(3) is an exact symmetry, and the 
constants values are taken from |19j. For 156 Gd the Huiq = 41 MeV xi -1 / 3 = T.QIQMeV is 
the oscillator constant for the equivalent spherical nucleus and the constants n and \x have the 
values = 0.06370,^ = 0.60, ^ = 0.06370, /i„ = 0.42 

4. Results and conclusions 

We studied the contribution of each term in the total Hamiltonian eq. Q by varying the 
strength of each interaction and plotting the energy of the first few states of angular momentum 
K = 0,2,4,6,8, and 1. In Figure [T] we present the energy values of few states of total angular 



momentum J = 0,2, and 4 versus the strength of the proton and neutron paring. In the 
calculation, the strength of the proton and neutron paring is set the same and varied. The 
paring interaction is important to get the correct relative value of the states in the K = 
band, versus the ground band. It is also important to get the relative ordering of the band head 
of the K = 2 + band versus the ground and K = 0^ bands. 




Figure 1. The energy values (in 
MeV) of the few three states of total 
angular momentum J = 0, 2, 4, versus 
the strength of the proton and neutron 
paring in 156 Gd. 
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Figure 2. The energy values in MeV, in the first 
four bands in 164 Dy. The first one is the ground 
band, the second one is the K = 2 + band, the third 
and fourth are the K = 0^~, and K 



0g bands. 



In the left hand side are given the experimental 
values[20j, and in the right the calculated ones. 



After varying one interaction strength at a time and studying the effect of each interaction on 
the energy levels and B(E2) transitions, we fixed the values for pairing (G w = 21/ A, G u = 17 /A), 
as well as for the quadrupole-quadrupole interaction strength (x = 35A~ 5 ^ 3 )[16\. The details 
of the study will be given elsewhere. The interaction strengths of the rotor terms were varied 
to give a best fit to the the band heads of the first excited K n = + , K T = 2 + and = 1 + 
bands, as well as the moment of inertia of the g.s. band. Explicitly, the term proportional to Kj 
breaks the SU(3) degeneracy of the different K bands, the J 2 term represents a small correction 
to fine tune the moment of inertia, and the last term, C2, is introduced to distinguish between 
SU (3) irreps with A and \x both even from the others with one or both odd, hence fine tuning 
the energy of the first excited K n = 1 + state. 

With parameters fixed, the Hamiltonian is diagonalized. As an example, the energy values 
for the first four energy bands in 164 Dy are given in Figure [2j In the left hand side are given 
the experimental values [20]. and in the right the calculated ones. Within this framework, the 
splitting and mixing of the pseudo-S'f7(3) irreps are generated by the proton and neutron single 
particle terms (H^l u ) and the pairing interactions. This mixing plays an important role in the 
reproduction of the behavior of the low- lying collective states in deformed nuclei |18j . 

This approach was employed for a systematic investigation of the experimental energies of 
the J w = 2+ J n = 0~t 1 , and J w = 2"J" . states in the deformed nuclei from the Fn = 

and Fo = 1 multiplets of the rare earths, which are shown in Figj3] and Fig. |4j where the 
rather irregular behavior of the band-heads of the first excited bands is quite obvious. The 
calculated energies in the framework of pseudo-SU(3) model, with above described fit of the 
respective parameter strengths, are compared to the experimental values in Tables [3] and |4j We 
calculated also all the states with J < 8 within these three bands. The obtained results are in 
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Figure 3. The experimental energy values 



[20] of the J n = 2+, ,r 







and 



if =27 



states in a sequence of nuclei 



for which the numbers of valence protons and 
neutrons are the same, (Fq = 0). 



Figure 4. The experimental energy values 



[20] of the J n 



2+ r 



K n =0. 



+ > 



and 



2^_ 2+ states in a sequence of nuclei 

for which the difference in the numbers of 
valence proton and neutron pairs is two, 
(F = 1). 



Table 3. Energy values for the J 77 = 27 , first excited J n = , and J n = 2 , states 

bJ 1 ' k*=QJ K 7r =2J 

for seven nuclei in the Fq = multiplet. The experimental values [20] are given in parenthesis. 
The numbers of valence proton and neutron pairs are given in the second and third columns. 
The resulted fitted parameters of the rotor part of the Hamiltonian are given in the fourth to 
seventh columns. 
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E(2 1 ) 
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Th. (Exp.) 


Th. (Exp.) 


Th. (Exp.) 














10" 4 


[MeV] 


[MeV] 


[MeV] 


152 Nd 


5 


5 


0.000 


0.00 


0.0000 


2.57 


0.082(0.073) 


1.14(1.14) 


1.31(1.38) 


156g m 


6 


6 


0.000 


0.55 


0.0000 


2.59 


0.078(0.076) 


1.07(1.07) 


1.45(1.47) 


160Q d 


7 


7 


0.001 


0.15 


0.0034 


1.93 


0.085(0.075) 


1.33(1.33) 


0.99(0.82) 


164 Dy 


8 


8 


-0.001 


0.04 


0.0008 


0.65 


0.073(0.073) 


1.67(1.66) 


0.76(0.76) 


168g r 


9 


9 


-0.002 


0.02 


0.0008 


0.75 


0.089(0.080) 


1.21(1.22) 


0.81(0.82) 


172 Yb 


10 


10 


-0.001 


0.12 


0.0010 


0.31 


0.081(0.078) 


1.04(1.04) 


1.49(1.47) 


176 Hf 


11 


11 


-0.007 


0.30 


0.0060 


0.43 


0.178(0.088) 


1.15(1.15) 


1.26 



very good agreement with experiment |18^ 120] . The main reason for obtaining the position of 
each collective band with respect to each other, as well as of each level within the band is the 
specific content of the obtained SU(3) irreps into the collective states, which is related to their 
deformations. For nuclei from Table [3] in the middle of the shell, the ground and the 7 band 
belong to the same (A, /x) with A > fi. At the limits of the deformed region the ground band 
states have oblate deformation (A < /i) and the K 77 = 0^~ and the K w = 2f bands are mixed in 
the same SU(3) irrep [18]. The analysis for nuclei from Table [4] is under investigation. 

The correct description of collective properties of first excited K n = Oj, and K w = 2 + states is 
a result of representation mixing and state deformation induced by the model Hamiltonian. This 
study shows that pseudo-spin zero neutron and proton configurations with relatively few pseudo- 



Table 4. Same as in Table [3] for four nuclei with .Fo = 1- 
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[MeV] 


[MeV] 


[MeV] 


156 Gd 


7 
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0.001 


0.10 


0.0014 


1.36 


0.094 (0.089) 


1.06 (1.05) 


1.14 (1.15) 


160 Dy 


8 
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0.003 


0.16 


0.0004 


0.94 


0.092 (0.087) 


1.30 (1.28) 


0.99 (0.97) 


168 Yb 


10 
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0.00 


0.01 


0.0010 


0.51 


0.088 (0.088) 


1.16 (1.15) 


0.98 (0.98) 


172 Hf 


11 
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0.00 


0.00 


0.0000 


0.51 


0.119 (0.095) 


0.87 (0.87) 


1.08 (1.07) 



SU(3) irreps with largest C2 values suffices to yield good agreement with known experimental 
energies. 

In the presented applications of the pseudo-SU(3) model, we achieved a good description of 
the complex behavior of a series of heavy deformed nuclei using a small configuration space. 
We are working on further analyzing the effect of the different components of the interaction 
in the Hamiltonian, on the strength of the transition probabilities. Also, the size of the model 
space will be further investigated. Work in progress is on expanding the study to other groups 
of nuclei related to one another by an over-arching F-type symmetry. 
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